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Abstract 

The role of electron correlation on different pairing symmetries are discussed 
in details where the electron correlation has been treated within the slave bo- 
son formalism. It is shown that for a pure s or pure d wave pairing symmetry, 
the electronic correlation suppresses the s wave gap magnitude (as well as the 
Tc) at a faster rate than that for the d wave gap. On the otherhand, a complex 
order parameter of the form (s + id) shows anomalous temperature depen- 
dence. For example, if the temperature (T^) at which the d wave component 
of the complex order parameter vanishes happens to be larger than that for 
the s wave component (T^) then the growth of the d wave component is ar- 
rested with the onset of the s wave component of the order parameter. In this 
mixed phase however, we find that the suppression in different components 
of the gap as well as the corresponding due to coulomb correlation are 
very sensitive to the relative pairing strengths of s and d channels as well as 



the underlying lattice. Interestingly enough, in such a scenario (for a case of 
Tjf > T^) the gap magnitude of the d wave component increases with electron 
correlation but not for certain values of electron correlation. However, this 
never happens in case of the s wave component. We also calculate the temper- 
ature dependence of the superconducting gap along both the high symmetry 
directions (F - M and F - X) in a mixed s + id symmetry pairing state and 
the thermal variation of the gap anisotropy ( Ar~ x ^ with electron correlation. 
The results are discussed with reference to experimental observations. 
PACS. No. : 74.25.DW, 74.62.-c, 74.20.Fg. 

I. INTRODUCTION 

The correlated electron systems like the high temperature cuprate superconductors, show 
various anomalous physical properties in the normal as well as in the superconducting state. 
Alongwith the unsolved problems concerning the pairing mechanism in these exotic mate- 
rials the question of order parameter symmetry is also not yet understood. In the weak 
coupling conventional BCS superconductors, superconductivity results from a pairing of 
electrons via phonon- mediated attractive electron-electron interaction which dominates the 
usual Coulomb repulsion at low temperature, and an energy gap A(A;) appears in the quasi- 
particle spectrum. This energy gap is very nearly isotropic in /c-space, so that the gap 
has the same magnitude and phase in all directions i.e., an isotropic s-wave pairing state. 
In the high-Tc superconductors the pairing of charge carriers is established by experiments 
on flux quantization I-V characteristics of Josephson tunnel junction p| and Andreev 
scattering 0. But the point contact tunneling on various high-Tc materials and nuclear 
quadrupole/magnetic resonance and relaxation measurements conclude that the energy 
gap function in these materials has considerable anisotropy. In high-Tc superconductors, 
measurements of the temperature dependence of the NMR Knight shift of suitable nu- 
clei in a superconductor ruled out the possibihty of p-wave pairing but whether the pairing 
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is of s-wave or d-wave type is not yet settled. The results of NMR experiments measuring 
Cu- relaxation rates by Martindale et al are in agreement with the prediction of d-wave 
pairing. Measurement of penetration depth A at low temperatures on YBa2Cu30Q_g5 sup- 
ports d-wave pairing. Recent angle - resolved photoemission spectroscopy (ARPES) study 
also suggests a d-wave state. However, there exist some experimental results like varia- 
tion of penetration depth in Nd2-xCexCuOi 0, measurements of Josephson supercurrent 
for tunneling between Pb and Y Ba2CuzO-j |]lO| which does not correspond to the exact 
d-wave symmetry. Moreover, some photoemission studies on Bi2Sr2CaCu20s+x [0] 



inconsistent with pure d-wave but are more consistent with a mixed state of s and d wave 
components. Kotliar [|12| and Ruckenstein et al, first introduced the concept of mixed 



s + id and s + d symmetries respectively. Their ideas of mixed configuration was to interpret 
the NMR and NQR data in the superconducting state of YBCO fl^ and the Joseph- 



son critical current measurements in YBCO (superconductor-normal-superconductor) SNS 
junctions and YBCO/Pb junction |T5[ . 



Based on the discussions on the experimental situation about the determination of the 
order parameter symmetry in high temperature superconductors it is clear that there are 
various opinions about pairing symmetry in cuprates, which are known to be correlated 
system. Therefore, in the present paper the effect of Coulomb correlation on the supercon- 
ducting state is considered based on a weak coupling theory. The role of electron correlation 
on different pairing symmetries are discussed in details where the electron correlation has 
been treated within the slave boson formalism. It is shown that for a pure s or pure d wave 
pairing symmetry, the electronic correlation suppresses the s wave gap magnitude (as well 
as the Tc) at a faster rate than that for the d wave gap. On the otherhand, a complex order 
parameter of the form [s + id) shows anomalous temperature dependence. For example, 
if the temperature {T^) at which the d wave component of the order parameter vanishes 
happens to be larger than that for the s wave component (Tj*) then the growth of the d 
wave component is arrested with the onset of the s wave component of the order parameter. 
In this mixed phase however, it is shown that the suppression in different components of the 



gap as well as the corresponding Tc are very sensitive to the relative pairing strengths of s 
and d channels as well as the underlying lattice. Interestingly enough, in such a scenario (for 
a case of > T^) the zero temperature gap magnitude of the d wave component increases 
with electron correlation but not for certain values of electron correlation. However, this 
never happens in case of the s wave component. We also calculate the temperature depen- 
dence of the superconducting gap along both the high symmetry directions (F - M and F - 
X) in a mixed s + id symmetry pairing state and the thermal variation of the gap anisotropy 
^ Ar~x ^ ^^^^ electron correlation. The results are discussed with reference to experimental 
observations. 

The lay out of the rest of the paper is as follows. In section II we describe the effect of on- 
site coulomb correlation on weak coupling superconductivity. In this section we also present 
the slave boson approximation to the on-site coulomb correlation in the Hubbard model. The 
next section HI is devoted to describe the effect of electronic correlation on different order 
parameter symmetries. Effect of electron correlation on each order parameter symmetry 
is described sequentially in each subsection of this section with detailed discussions on the 
results. A comparison of the present study with ARPES results in high-Tc materials are 
described in section IV. Finally, remarks have been made about the present calculation in 
the conclusion section V. 



II. EFFECT OF ELECTRON CORRELATION ON A WEAK COUPLING 
THEORY OF SUPERCONDUCTIVITY 

The simplest model that describes the electron-electron interaction in a correlated sys- 
tem, like the high-Tc cuprates, is the single band Hubbard model. One of the major problems 
in solving the Hubbard model is how to treat the correlation exactly. One of the ways to 
do so is to take recourse to the slave boson formulation of Kotliar and Ruckenstein (KR) 
p!8| . In this representation a set of four bosons are assigned for four possible occupancies of 
a lattice site. The bosonic fields which keep track of the four different occupations of a site 
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'i' are e|(ei), sJ^/(_^-)(sjo-/(_o-)) and d\{di) corresponding to creation (annihilation) of empty, 
single occupation with spin cr/(— a) and double occupation respectively. Of course the total 
probability of occupation of any site is one and it has to be respected. Secondly, the fermion 
number should be conserved at 'i'th site in this slave boson representation. These two 
constraints are imposed by two conditions, completeness relation and charge conservation, 
written as, 

4ei + J2 -^t-Sia + 4di = 1 (1) 

^icT^i(T — ^\a^io- + d\di (2) 

The Hubbard Hamiltonian is described as 

ijcr % id 

where the first term is the usual tight binding Hamiltonian with the hopping integral tjj, the 
second term denotes the intra-atomic Coulomb repulsion (f/) between two electrons with 
opposite spins occupying the same site and the third term denotes the chemical potential 
(yu) which can take into account the deviation from half - filling. The prescription for the 
transformation of equation (|^) in terms of the slave boson operators is given below : 



(4) 



where 

= (1 - dldi - sl^Si^)~^{elsia + sl^^di){l - efcj - sj-a-Si-a)"^ (5) 

and cfcr is the modified fermion creation operator. The form of Zi„ is so chosen as to reproduce 
the correct band structure in the absence of correlation. The operators (efsjo- + sl_^di) in 
the eqn. (^) describe the hopping process of electron i.e., if an electron hops from site 'i' to 
'j', the slave bosons must simultaneously change at j and i. Depending on whether the site 



'i' is singly or doubly occupied the bosonic state of 'i' must change from to e| or from 
dl to sj_^. Thus there are two transition channels which add up, and the total transition 
probability must be equal to one. It is therefore useful to introduce a normalization factor, 
which guarantees the conservation of probability even in the mean-field theory. 

Following KR the Hubbard Hamiltonian can be re-written in terms of the slave boson 
operators as 



ija i icr 

icr 

+ Yl Ki^ ~ ~ ~ Yj ^\<jSic) (6) 

i cr 

where Ajo- and are Lagrange multipliers which enforce the charge conservation and com- 
pleteness relation (1,2) respectively. So the Coulomb interaction term is simplified and 
diagonalized with double occupancy operators but the kinetic energy part becomes compli- 
cated. The values of the boson field operators and the Lagrange multipliers are determined 
by minimizing the free energy of the system in the saddle point approximation, where all the 
bose operators and Lagrange multipliers are replaced by c-numbers. In this approximation 
the Hamiltonian is 



H = J2 QtijclaCja + {X- ^iaCia 
ija ia 

+N \Ud^ - 2X{d^ + s^) + A'(l -d^-e^- 2s^)] (7) 

where q =< zj^Zja- >, e, d, s, A, A' arc the saddle point values of the respective field operators, 
W is the number of sites, and in the case of paramagnetic ground state Sa- — S-o- = s. The 
Hamiltonian of the system takes the simple form of an effective tight binding model with 
a modified hopping integral of the form tij — > qtij, where the correlation effects arc built 
in through the multiplicative factor q —< z^z >. q in general is a complicated function 
of the coulomb correlation u{— Uc being the Brinkmann-Rice value for metal-insulator 
transition) and the dopant concentration 6. In this approach solutions are obtained for 



the paramagnetic states for all values of u and band fillings that reproduces the correct 
Brinkmann-Rice result for metal-insulator transition at a critical value of correlation (?7c) 
at half-filling. Approximately, in the weak correlation limit q = l — u^ and in case of strong 
correlation and small values of 5, g = ^^^^^^ (for details please see ref. 

Now in order to study the effect of electron correlation on SC pairing symmetry we use a 
model Hamiltonian which contains in addition to a repulsive on-site coulomb correlation (7) 
term, a pairing term that leads to superconductivity (for our purpose). The pairing interac- 
tion is assumed to be due to some boson exchange (unknown) mechanism, as there exists no 
conclusive pairing mechanism for the high-T^ superconductivity. The pairing hamiltonian 
in general may be written as, 

Hp= -Y. Vkk'cl^cl^^c.k'iCk"\ (8) 

kk' 

While treating the electron correlations using salve boson formalism, the nature of the 
quasiparticles are no longer the bare electrons but the electronic quasiparticles of the corre- 
lated system (i.e c). Therefore the bare electronic operators should be transformed accord- 
ingly as equation (4). As a consequence, the strength of the pairing interaction becomes 
Vkk' — ^ q'^Vkk' when the Cka{cl„) are replaced by new quasiparticle operators Cka{cl„). 
Hence, the total hamiltonian for the superconducting state in a correlated system can be 
written as, 

H = Ho + Hn + Hp (9) 

where 

^0 = I](gefe - /i)cLcfea (10) 

ka 

Hb = UY^ d\di - ^ \ia{s\„Sia + d\di) + ^ A-(l - ejcj - ^ s\„Si„ - d\di) (11) 

i ia i o- 

Hp = -Y. ^kk'q^c{^c^-k[C-k'iCk'] (12) 

kk' 



19,20|). 
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where, = —2t [{cosk^a + '^icoskyO) + 'y2Coskj.acoskya\, where 72 = t and t' represent 
nearest neighbour and next nearest neighbour hopping respectively, 71 = 1 for square lattice 
and in presence of orthorhombic distortion 71 < 1. 

Since in the saddle point approximation the bosonic part Hb is constant, so the effective 
Hamiltonian is given by 

H = Ho + Hp (13) 

It can be seen from eqn. ([T2|) that the pairing vertex between the correlated electrons 
mediated by some bosonic exchange is renormalized to q'^Vkk'- Since deviates from unity 
(g^ < 1 for M 7^ 0), the SC-pairing amplitude will always be suppressed in presence of 
electron correlation, however, its details will depend on the value of u, 6 and the nature of 
Vfc fc/, is the main point of investigation in this work. 

The SC order parameter in a correlated system may be defined as 

A{k) = Y.q'Vkk'<5l^5l,,^> (14) 

k' 

The Hamiltonian ( |13|) is treated within the mean field theory (|I^) in order to yield the SC 



gap equation within a weak coupling theory as. 



Aik) = Y.fV,,,^t^nhA (15) 

where the quasiparticle energy is given by 

El = iqe,-^^r+\Aik)\' (16) 

In the next section we will discuss about the different pairing symmetry and its effect in a 
correlated system. 

III. EFFECT OF ELECTRON CORRELATION ON ORDER PARAMETER 

SYMMETRY 



Based on the nature of pairing potential the symmetry of SC gap could be different. The 
pairing potential is usually assumed to have a separable form for simplicity i.e., 
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Vkk' = V7]kT]k' 



(17) 



Depending on the nature of the /c-dependence of r^^ one gets different symmetries such as 
(i) r]k = constant, corresponds to an isotropic s-wave (conventional BCS) symmetry, (ii) 
Vk = f{k), refers to an anisotropic s-wave symmetry, if f{k) is a smooth function in the first 
Brillouin zone and is positive definite (i.e., nodeless), (iii) r^^ = cosk^a + coskya, corresponds 
to an extended s-wave pairing symmetry and (iv) r/^ = coskxa — coskya corresponds to 
dx2_y2 pairing symmetry. In the following we shall discuss only about the isotropic s-wave, 
the (i-wave and a mixed symmetry state which will be defined later on. 



and the corresponding SC gap function = A^. Hence the gap equation for a pure s-wave 



The gap equation (\[% can be solved analytically following the standard procedure of the 
BCS gap equation. However, we solve the gap equation numerically as the same will not be 
possible to solve analytically in case of (i-wave as well as for the complex order parameter 
(discussed later on). All the parameters are expressed in units of 't'. However, at the end 
temperature has been expressed in Kelvin assuming t = 0.3 eV (realistic for copper oxide 
system), in order to examine, the real transition temperatures. For a set of parameter 
values Uc = 0.6 (cut-off value), Vg = 1.0 we present the numerical results in Fig. 1 and Fig. 
2. Fig. 1 describes the temperature variation of the pure s-wave gap for a perfect square 
lattice whereas the Fig. 2 represents that for a distorted lattice (tx ^ ty , i.e., 71 = 0.95) 
with next nearest neighbouring hopping t' = 0.02t. It is evident from both the figures 1 
and 2 that the magnitude of A^ as well as the corresponding Tc decreases sharply with u. 



A. Pure s-wave 



The form of the pairing potential for isotropic s-wave symmetry is given by 



Vkk' = Vs = constant 



(18) 




(19) 
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Furthermore, for the same set of parameter not only that the magnitude of and % is less 
in the distorted case (cf. Fig. 2) than that for the perfect square lattice but also the fall in 
Tc with increasing u is more in comparison to the former case. So, the suppressive effect of 
electronic correlation on and Tc is comparatively more in the orthorhombic phase (this 
may be clearer from Fig. 5) as will be discussed in the later subsections. 

B. Pure d-wave 

In case of the pure d^2_y2 symmetry the pairing potential is given by 

Vkk' — Vdicoskx — P'cosky){cosk'^ — f3'cosk'y) (20) 

and the corresponding gap function is assumed as, 

^d{k) = A°(cosA;, - (5'cosky) (21) 

where Vd and A^J represents the strength of the li-wave pairing potential and the amplitude 
of the gap function respectively. Where (3' — 1 or < 1 corresponding to a square lattice or 
a orthorhombic lattice. So the resulting gap equation is written as 

„ „ ^^PJcosk'- l3'cosk'f ^ ,3Ek', , , 

= fva E — ^rr^^ tanh(^) (22) 

In Fig. 3 and Fig. 4 we show the thermal variation of the d-wave gap for the set of parameter 
values ojc — 0.6, — 0.5 for square (71 = /3' = 1) and orthorhombic (71 — (3' — 0.95) 
lattices with next nearest neighbour hopping t' = 0.02t respectively. The qualitative trend 
of behaviors for the thermal variation of the gap parameter remains almost unaltered i.e, 
the magnitude of the gap as well as the Tc reduces with u. In case of orthorhombic lattice 
the SC gap and Tc is smaller than that of the undistorted one. 

Fig. 5 shows the variation of T° — Tc vs u for s and d wave symmetry in distorted and 
undistorted lattice. T° is the SC Tc for u — 0. It is evident from Fig. 5 that for distorted 
lattice the effect of u on Tc is more than the undistorted one for both s and d wave symmetry. 
It may be compared from figures 1 and 3 that for an increase in u from to 0.4 there is a 
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reduction in Tc about 50% and 44% for the s and the d wave gaps respectively. The same 
when compared between the figures 2 and 4 the reduction in is about 58% and 54% for the 
s and (i-wave respectively. Therefore, in general the coulomb correlation affects the s-wave 
symmetry largely in comparison to the d-wave gap. 



C. Mixed s + id order parameter symmetry 

In the mixed symmetry state there is a mixture of both s and d wave order parameter and 
a phase difference of e*^ is considered The generalized form of the pairing interaction 
causing the superconductivity is given by 

Vkk' = Vs + Vd{coskx — P'cosky){cosk'^ — (3'cosk'y) (23) 

where Vs and Vd corresponds to the strengths of s and d wave channel interactions. And the 
corresponding s + id wave order parameter is taken as, 

A(A;) = + iA2(cosfc^ - (3'cosky) (24) 

Substituting the equation (pSf) and (|2^) into equation (|15D and separating the real and 



imaginary parts of the equation we obtain gap equations for coupled s and d wave compo- 
nents given by 

A, = q^v, V tanh(^^^) (25) 

k' ^^k' ^ 

and 

\ = <i^dl^ ^t;— ^tanh(^— ) (26) 

with 

i?r"' = {q^k - /i)' + A^ + Af{cosk, - (3'coskyf (27) 

Self-consistent solution of equation ( [25| ) and ( PB| ) gives A^ and A° in a mixed (s + id) wave 
state and for either t>d = or f ^ = the coupled gap equations reduce to pure s-wave or pure 
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(i-wave (as discussed in the earlier subsections) respectively. To note, the difference between 
(p!9|) and ( |25D as well as that between (p2D and (^) lies mainly with the difference between 



Ek and E'^'^^. Here we performed the self-consistent solution of ( |25D and ([26| ) with a cut 
off LOc = 0.6 (as in the earlier cases of a pure s or a d wave scenario). 

In case of mixed s + id symmetry with Vd and Vg both being finite, Fig.6 shows the 
phase diagram where the magnitudes of the order parameters with different symmetry at 
a temperature of lOK are plotted against the ratio of the strength of interactions (^) for 
different u for the square lattice. This phase diagram has three regions, the region where 
only s-wave solution exists i.e, for ^ < 0.45, only d-wave solution exists i.e, for ^ > 0.39 and 
both the s and d wave solution co-exist for 0.39 < — < 0.45. For finite electron correlation 
u = 0.4 the region with mixed symmetry shrinks to within 0.39 < — < 0.42 with larger 
decrease in the magnitude of the s wave gap compared to that for the rf-wave component. 
To note, with the electron correlation increased (to u = 0.4) the region (in terms of ^) 
where the s -wave solution exists reduces (to ^ < 0.42), but that for the rf-wave remains 
unchanged although the magnitude is suppressed. It is also important to note that the 
magnitudes of the s or d- wave component gaps with correlation depends crucially on the 
relative strength of interaction — . 

Furthermore for ^ ~ 0.39 the s-wave gap decreases with u while ci-wave gap increases 
slightly from its value compared to that at the uncorrelated case. This is because the 
rate of decrease of s-wave gap is more rapid than the d wave gap with correlation. For a 
clearer observation, the variation in the temperature dependence of the s-wave and d-wave 
component of the gap parameters in the mixed state (for ^ = 0.4) with electronic correlation 
is shown in Fig. 7. The graph with diamond shows pure s-wave gap for the uncorrelated 
case [vg = 1, = 0). In the mixed state ^ = 0.4, (and in the absence of correlation) 
the s-wave gap has value slightly less than that of the pure case (the dashed line). With 
increasing correlation, u = 0.4, while the magnitude of s-wave gap is suppressed strongly 
(dotted line), that of the d-wave gap acquires higher value than the u = case (the solid 
line). Nevertheless, the d-wave (TJ) is higher for the uncorrelated case than that with 
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correlation. Hence, in the lower temperature regions there is a competing effect between the 
s and c?-wave gap parameters in the mixed symmetry region; with the appearance of the 
s-wave order parameter the d-wave gap is suppressed for the u = 0.4 case (the dash dotted 
line) . The d-wecve gap attains its highest value where the s-wave component gap amplitude 
vanishes. In this figure, the curve with +-1- sign represents the pure c?-wave gap with all the 
parameters being same {vd = 0.4) ; demonstrating the arresting of the growth of the d-wave 
component with the opening up of the s-component gap in the mixed phase. 

Figures 8 and 9 show the variation of s-wave and d-wave gap parameters for different 
u as the ratio of the strengths of interactions increases close to and beyond the values 
corresponding to the co-existence regions (for u — 0.4) i.e., at ^ = 0.42 and ^ = 0.45 
respectively. The gap magnitude for both the order parameter symmetry decreases with 
electron correlation, however, the rate of decrease is different for s and d in the mixed state 
which depends on ^ ratio as well as the electron correlation. Here for the uncorrelated case 
it is clearly visible that the growth of d-wecve component of the gaps is arrested with the 
onset of the s-wave component gap. In fact a detailed investigations of the Figures 7 to 9 will 
demonstrate that in the mixed phase within the coexistence regime i.e (0.39 < ^ < 0.42), 
the suppression in different components of the gap as well as the corresponding T^s are very 
sensitive to the relative pairing strengths of s and d channels. 

So far all the calculations for the mixed phase, were carried out for the square lat- 
tice. However, it was seen in case of the pure s and d phases that the introduction of 
orthorhombicity causes reduction in both, the magnitude of the order parameters as well 
as the corresponding TcS. So the effect of orthorhombicity on the mixed phase is now in- 
vestigated. The phase diagram for different order parameter symmetry as a function of the 
^ ratio is depicted in the Fig. 10 for the orthorhombic distortion and finite next nearest 
neighbouring interaction. The s-wave and d-wecve gap has lesser magnitude in the distorted 
lattice than that for the perfect square lattice in mixed s + id state and the region of s and d 
co-existence is shifted to smaller — value and is shrunk to smaller region of area compared 
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to the undistorted lattice. Similar to the phase diagram in Fig. 6, the value of ^ for the 
occurence of d-wave component does not change with u whereas that for the s component 
changes from ^ < 0.48 (for m = 0) to ^ < 0.46 (for u = 0.4). With a higher value of 
Coulomb correlation {u — 0.4) the s — d mixing region reduces further along with a large 
suppression in the magnitude of s-wave and d-wave gap parameters. It is noticable that 
the s-wave component suffers larger suppression in magnitude of the gap than that for the 
d-wave. The on-site Coulomb correlation strength u as well as the isotropic pairing strength 
Vs both being isotropic in nature, the s-wave gap magnitude is largely affected by Coulomb 
correlation, whereas due to extended nature of the d-wave component it is less affected by 
correlation. However, we note that in the mixed phase, depending on the value of ^ the 
effect of electron correlation on transition temperatures or is very different. This is 
due to the fact that in the mixed phase none of the components (s or d) follow BCS like 
temperature variation. This has sometime led to increase in the gap magnitude with u but 
reduction in the (= T^) (cf. Fig. 7). 

Fig. 11 shows the variation of s and d wave component with electron correlations for 
^ = 0.45. Similar to the square lattice the s and d wave gap parameters in the pure state 
take higher values than that in the mixed state for same parameter values v^,, Vg and u. 
This figure has to be contrasted with Fig. 9. With the introduction of orthorhombicity the 
nature of the temperature variations of the different gap components changes abruptly. It is 
worth noting that in the pure s or pure d wave cases (discussed earlier) the orthorhombicity 
suppresses both the magnitude of the order parameters and the corresponding T^s (cf. Fig. 1 
to Fig. 4) . However in the mixed state, because of the interplay between the two components 
of the order parameters, the orthorhombicity together with next nearest neighbour hopping 
enhances the relative magnitude of the s-wave component whereas reduces that of the ci-wave 
when compared with respect to the perfect square lattice for same value of ^ (cf. figures 9 
and 11). Considering the case of it = 0.4 from Fig. 9 and Fig. 11 it would be self-evident that 
the in the orthorhombic phase reduces very largely from its value for the square lattice, 
whereas in the orthorhombic phase increases strongly. On the other hand unlike the 
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undistorted lattice, with the mixed symmetry (cf. Fig. 6), where the d-wave component can 
have higher magnitude for higher u close to ^=0.4 ; in presence of orthorhombic distortion 
the (i-wave order parameter never acquire higher magnitude for higher u, for any value of 
^, as is clear from Fig. 10 as well as Fig. 11. 

IV. TEMPERATURE DEPENDENT GAP ANISOTROPY IN A CORRELATED 

SYSTEM 

There are three broad catagories of experiments, which had been used to probe the 
symmetry of the superconducting gap function ||8|- [l0| , p!6l , ^^ namely (i) magnetization and 



transport, (ii) spectroscopies and excitations and (iii) Josephson measurements. The most 
direct and fundamental probe of the gap magnitude in superconductors is angle resolved 
photoemission spectroscopy (ARPES), which indicates a highly anisotropic gap in high Tc 
compounds. 

The main features observed from ARPES by different groups of experimentalists are 
summarized below : Shen et al, observed that the SC-gap has (a) maximum magnitude 
along r — M direction (Cu-0 bond direction in real space) ; (b) minimum magnitude (neg- 
ligible spectral weight) along P — X(Y) direction (diagonal to the Cu-0 bond) and there is 
a monotonic increase in the magnitude of the gap from its smallest value along the P — X 



direction, thereby confirming true anisotropy. The authors of ||23| argue that these results 



are consistent with d^'^-y'^ symmetry of the gap. However, the data on the momentum 



resolved temperature dependence of the SC gap of Bi-2212 by J. Ma et al [11] is not only 



in contradiction with Shen et al p3[ but also it shows some exciting features. For mea- 



surements of the temperature dependence of the gap, the angle resolved photoemission data 
have been taken for temperatures ranging from 36 to 95 K along the P — M as well as in the 
P — X (40 to 95 K) high symmetry directions. The photoemission SC condensate spectral 
area along P — M direction exists till the bulk transition temperature, Tc = 83K. The mea- 
sured SC gap is almost independent of the temperature at lower temperatures and retains 
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its maximum value even up to temperature 85% of its Tc. But most interesting features to 
note is that, the photoemission condensate spectral area is remarkably weaker at 40K in the 
r — X direction compared to that in the T — M direction, indicating a smaller size of the 
gap. The size of the gap at 36 K in the T — M direction is 16 meV whereas that in the 
r — X direction is 10 meV (at T=40 K). The finite gap magnitude along F — X direction 
below Tc ruled out the possibility of simple rf-wave pairing [|l6l. However, at 70 K when the 
gap along the T — M direction retains 90 to 100% of its full value, becomes indistinguishable 
from zero along the T — X line. This may be taken as a signature of a two component order 
parameter, dx2_y2 type close to Tc and a mixture of the both s and d otherwise. Moreover, 
the significant point is that, due to the unconventional temperature dependence of the SC 
gap in different parts of the Brilloun zone, the temperature dependent gap anisotropy i.e., 
the ratio ^^'^(^~-^) jg enhanced closed to Tc. In the present section we show that the above 
peculiarity in the temperature dependence of the SC gap is a natural consequences of the 
mixed s + id symmetry and in a correlated system the gap anisotropy increases with electron 
correlation. 

In the gap equations (E3) and (EBF) and A° represents the amplitudes of the s and d 



wave gaps respectively. Since represents the isotropic s-wave gap so it retains its constant 
maximum value for any kx and ky direction whereas magnitude of A^ varies with k^ and ky 
direction due to the multiplicative factor {cosk^ — cosky). In the F — M(M) direction i.e., 
(±7r, 0), (0,±7r), {coskx — cosky) acquires maximum magnitude =f2, while in the F — X{Y) 
direction (±f , ±f ) {cosk^ — cosky) is zero identically. So it is expected that the rf-wave gap 
parameter will attain its maximum at F — M{M) points and minimum at F — X(Y) points, 
however s-wave gap is constant all over the Fermi surface. Considering the energy spectrum 
of the mixed phase EJ!^^^ we write A(F — M) and A(F — X) as 



A(F-M) = VA2 + 4A2 (28) 

and 

A(F-X) = A, (29) 
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Hence the ratio which is a measure of gap anisotropy can be written as 



A(r - M) _ y/A^ + 4Ai 

A(r-x) - ^^"^ 

Fig. 12 shows the variation of A(r — M) and A(r — X) for u = and u = 0.4 in an 
orthorhombic lattice with ^ = 0.46. The A(r — M) is larger and almost constant with 
temperature till A(r — X) is zero, then A(r — M) falls to zero at the T^. The decrease of 
A(r — M) with u is much less than that for the A(r — X). The inset figure of Fig. 12 
depicts the temperature dependent gap anisotropy in a correlated case. The ratio ^^jZ.^) 
is almost constant at lower temperatures and increases faster near the temperature where 
A(r — X) = As is going to be zero. ratio is higher for higher correlation. This 

feature is however not an unique consequence of s + id pairing symmetry e.g, similar results 



are also obtained in a modified version of spin bag [24] mechanism of superconductivity. 



V. CONCLUSION 

We have presented a detailed study of the effect of electron correlation on pairing symme- 
try within a weak coupling theory. The method used to treat electron correlation is the slave 
boson formalism of KR [18] which reproduces the Brinkmann-Rice value of metal-insulator 
transition (Uc) correctly in the paramagnetic state. The value of u used in this work is 
always less than 1 (i.e U < Uc) so that weak coupling mean field theory of superconductivity 
can be applied. All through out the present calculation we restricted to half-filled situa- 
tion. Within this study we found that electron correlation has important effect on pairing 
symmetry. The detailed nature of order parameter, its temperature variation as well as 
magnitude depends very sensitively with electron correlation and the nature of underlying 
lattice considered. Few interesting features of interplay between the different components of 
the order parameter with electron correlation has been demonstrated (in the s + id picture). 
Such study may have some bearings to the high temperature superconductors. However, 
we have made no effort to understand the mechanism of superconductivity for cuprates in 
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this work. The electron correlation always affect the pure s wave gap more than the pure 
d wave symmetry gap irrespective of the underlying lattice. In the mixed s + id phase, the 
suppression in different components of the gap as well as the corresponding Tc are very sen- 
sitive to ^ and the underlying lattice. Assuming the s + id nature of the order parameter, 
the measured temperature dependent gap anisotropy from the ARPES of J. Ma et al, [11] 
can be understood on which the role of electron correlation has been emphasised. However, 
such a picture is not unique [24] . 
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Figure Captions 



Fig. 1 Thermal variation of the isotropic s-wave gap {vs = 1) for various values of the 
electron correlation {u — ^) in a perfect square lattice (for other parameters, see text). 

Fig. 2 Thermal variation of the isotropic s-wave gap {vg = 1) in an orthorhombic lattice 
for various values of u (all other parameters being same as that in Fig. 1). 

Fig. 3 Thermal variation of the pure d-wave gap {vd = 0.5) for various values of the electron 
correlation in a perfect square lattice (for other parameters see text) . 

Fig. 4 Thermal variation of the pure d-wave gap {vd = 0.5) for various values of electron 
correlation in an orthorhombic lattice (with all the parameters being same as Fig. 3). 

Fig. 5 Change in Tc due to correlation for pure-s and pure-d wave gaps (in square and 
orthorhombic lattices) where = Tc{u — 0). 

Fig. 6 Phase diagram of a s + io? superconductor, where the amplitudes of the s-wave and 
d-wave component of the order parameter in a square lattice is plotted as a function of their 
relative strength ^ for different values of electron correlation at T = 10 K. 

Fig. 7 Thermal variation of the amplitudes and A^; (in units of t) of the complex order 
parameter (s -|- id) for different values of electron correlation at ^ = 0.4 in a perfect square 
lattice (all other parameters are kept fixed as earlier figures) . 

Fig. 8 Thermal variation of A^ and A^ of the complex order parameter (s + id) for different 

values of electron correlation at ^ = 0.42 in a perfect square lattice (all other parameters 
are same as earlier figures). 
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Fig. 9 Thermal variation of and of the complex order parameter (s + id) for different 
values of electron correlation at ^ = 0.45 in a perfect square lattice. 

Fig. 10 Phase diagram of a s + id superconductor where the amplitudes s-wave and d-wecve 
component of the order parameter is drawn against their relative pairing strength in an 
orthorhombic lattice for different electron correlations (T = lOK). 

Fig. 1 1 Thermal variation of A^ and A^ of the complex order parameter (s + id) for differ- 
ent values of electron correlation in an orthorhombic lattice at ^ = 0.45. 

Fig. 12 Temperature variation of the superconducting gap along the two high symmetry 
directions for u — 0.0, 0.4 in the mixed s + id phase. The inset figure depicts the ratio of the 
superconducting gap along the two high symmetry directions ( a"" y-xj ) '^ith temperature 
for u = 0.0, 0.4. 
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